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POPULAR SUMMARY 

Satellite measurements at wavelengths of light corresponding to about 3.7 
jum have been exploited for almost three decades. Radiation near this wavelength 
(which includes both solar energy and energy emitted by the Earth's surface and 
atmosphere) occurs in a so-called "atmospheric window" where the atmosphere 

is relatively transparent to such light. Applications are primarily related to cloud 

< / ■ 

studies, including doud detection, doud thermodynamic phase, and quantitative 
doud particle size retrievals. The utility of the band has led to the incorporation 
of such channels on a number of existing satellite imagers (AVHRR instrument 
on the NOAA polar orbiters, MODIS instrument on NASA Terra and Aqua 
platforms, as well as others) and future operational imagers (VIIRS instrument 
on NPOESS and the precursor NASA NPOESS Preparatory Project (NPP) 
mission). 

Fundamentally, cloud information is derived by first converting the 
measured signal into a reflectance quantity. Conversion of calibrated 3.7 jum 
channel measurements (in radiance units) to reflectance requires knowledge of 
the incoming solar energy (spedfically the irradiance) in this spectral region. 
Despite the ubiquity of 3.7 pm channels, absolute solar spectral irradiance data 
comes from either a single measurement campaign in 1969 or models that predict 



the solar irradiance spectrum. 

In this study, we compare historical 3,7 pm band spectral irradiance data sets 
with the new semi-empirical solar model of the quiet-Sun by Fontenla et al. 
(2006). The model has expected uncertainties of about 2% in the 3.7 gm spectral 
region. We find that 3.7 jum channel-averaged spectral irradiances using the 
Thekaekara et al. observations or previous solar models differ from the new 
model by -1.5 to +4.1% (depending on the particular AVHRR and MODIS 
instruments). For MODIS, these solar irradiance uncertainties result in doud 
microphysical retrievals uncertainties comparable with other fundamental 
reflectance error sources. 
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ABSTRACT 

Since the launch of the first Advanced Very High Resolution Radiometer 
(AVHRR) instrument aboard TIROS-N, measurements in the 3.7 |rm atmospheric 
window have been exploited for use in cloud detection and screening, cloud 
thermodynamic phase and surface snow/ice discrimination, and quantitative 
cloud particle size retrievals. The utility of the band has led to the incorporation 
of similar channels on a number of existing satellite imagers and future 
operational imagers. Daytime observations in the band include both reflected 
solar and thermal emission energy. Since 3.7 pm channels are calibrated to a 
radiance scale (via onboard blackbodies), knowledge of the top-of-atmosphere 
solar irradiance in the spectral region is required to infer reflectance. Despite the 
ubiquity of 3.7 pm channels, absolute solar spectral irradiance data comes from 
either a single measurement campaign (Thekaekara et al. 1969) or synthetic 
spectra. 

In this study, we compare historical 3.7 pm band spectral irradiance data sets 
with the recent semi-empirical solar model of the quiet-Sun by Fontenla et al. 
(2006). The model has expected uncertainties of about 2% in the 3.7 pm spectral 
region. We find that channel-averaged spectral irradiances using the 
observations reported by Thekaekara et al. are 3.2-4.1% greater than those 
derived from the Fontenla et al. model for MODIS and AVHRR instrument 
bandpasses; the Kurucz spectrum (1995), as included in the MODTRAN4 
distribution, gives channel-averaged irradiances 1.2-1. 5 % smaller than the 
Fontenla model. For the MODIS instrument, these solar irradiance uncertainties 
result in cloud microphysical retrievals uncertainties comparable with other 
fundamental reflectance error sources. 
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I. Introduction 

The 3.7 pm atmospheric window was first used in satellite Earth remote 
sensing with the 1978 launch of the original 4-channel Advanced Very High 
Resolution Radiometer (AVHRR) instrument aboard the TIROG-N polar orbiter. 
AVHRR has retained this spectral capability up to the present (currently through 
the NOAA-18 polar platform, though post NOAA-15 platforms carry the 
AVHRR/ 3 version of the instrument which has a single data channel shared by 
both 3.7 pm and 1.6 pm spectral channels). The main objective in originally flying 
the channel was reported to be for cloud screening of sea surface temperature 
observations (Schwalb 1982). Measurements in this window are now routinely 
used for cloud detection and screening (Saunders and Kriebel 1988; Ackerman et 
al. 1998; Heidinger et al. 2002); fire detection (Kaufman 1991; Prins and Menzel 
1992, Justice et al. 2002); cloud phase and surface snow/ ice discrimination 
(Pavolonis et al. 2005); and quantitative cloud microphysical retrievals (Arking 
and Childs 1985; Platnick and Twomey 1994; Han et al. 1994; Minnis et al. 1995; 
Platnick et al. 2003). The usefulness of the band for cloud observations essentially 
derives from the significant dependence of single scattering albedo on cloud 
thermodynamic phase and particle size (absorption increases for the ice phase 
and with particles size), and differences in single scattering albedo (and thereby 
cloud emissivity) compared to window IR channels. Smaller land surface 
emissivity compared to window bands also provides contrast. Qualitative 
overviews of 3.7 pm imagery from AVHRR and its uses for clouds and sea ice 
discrimination can be found in Scorer (1986, 1989). 

The usefulness of the spectral band has led to the incorporation of a 3.7 pm 
channel on a number of low Earth orbit imagers over the last decade. Examples 
include the Along-Track Scanning Radiometer (ATSR) on board ERS-1/-2 and 
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EVNISAT, Visible Infrared Scanner (VIRS) on the TRMM spacecraft (Barnes et al. 
2000), Moderate Resolution Imaging Spectroradiometer (MODIS) on the Terra 
and Aqua spacecrafts (Barnes et al. 1998), and Global Imager (GLI) on Midori-II 
(ADEOS-II) (Nakajima et al. 1998). On geosynchronous platforms, 3.9 pm 
channels are available on the new GOES Imager and the Meteosat Second 
Generation (MSG) Spinning Enhanced Visible and Infrared (SEVIRI) instrument 
(Aminou 2002). Future operational imagers will continue to fly similar channels 
for weather and climate applications, including the Visible Infrared Imager 
Radiometer Suite (VIIRS) on the National Polar-orbiting Operational Satellite 
System (NPOESS) and NASA's NPOESS Preparatory Project (NPP) (Lee et al. 
2006), AVHRR/3 on the European Space Agency MetOp platforms, and the 
GOES-R Advanced Baseline Imager (ABI) (Schmit et al. 2005). Similar channels 
have also been included on aircraft imagers (e.g., the MODIS Airborne Simulator 
flown on NASA's high altitude ER-2 aircraft. King et al. 1996). 

Daytime quantitative use of the 3.7 pm band for cloud microphysical 
retrievals must account for contributions from both reflected solar and thermal 
emission energy. The observed cloud bidirectional reflectance and emissivity, 
and not the measured radiance itself, are the fundamental quantities relevant for 
cloud microphysical retrievals. In this paper, we examine the fundamental 
uncertainty in the observed 3.7 pm solar reflectance that corresponds to 
uncertainty in the knowledge of the solar irradiance in this spectral region. The 
partitioning between solar and emissive radiance in practical cloud retrieval 
algorithms is not discussed here. 

While many modern imagers have onboard solar diffusers to provide a 
direct reflectance calibration scale for solar spectral channels (e.g., MODIS), these 
calibration systems do not extend into the 3.7 pm region. A 3.7 pm channel is 
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calibrated to a radiance scale (via onboard blackbodies) and knowledge of the 
top-of-atmosphere (TOA) solar irradiance across the channels spectral bandpass 
is required to infer reflectance. The observed bidirectional reflectance of an 
Earth-atmosphere scene is 


R = 


31. 1 

MqFq 


( 1 ) 


where no is the cosine of the solar zenith angle, I is the satellite measured 
radiance for the viewing geometry, and Fo is the TOA solar irradiance. All 
radiative quantities in Eq. (1) are averages over the spectral bandpass of the 
instrument, i.e.. 


/ 

0 

/,/(A)J»(A)cg l 

J a F 0 (Xy!>(X)dX 


( 2 ) 


where <£>(A) is the instrument channel spectral response function. The relative 
uncertainty in reflectance due to uncertainty in the averaged irradiance over the 
channel can be written as 
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with Foi and F02 representing calculations from two different spectral irradiance 
datasets. Note the different denominator indexes in the exact expression. As an 
example, an overestimation in the band-averaged solar irradiance results in an 
underestimation of reflectance. 

Despite the widespread use of 3.7 pm satellite and aircraft measurements, 
only a single set of observations of absolute TOA spectral irradiance from the 
entire solar disk has been reported in this spectral region (Thekaekara et al. 1969; 
Thekaekara 1974). Further, the measurements are at a relatively poor spectral 
resolution. Here we give an overview of the historical observational data and 
common synthetic data sets of spectral irradiance in the 3.7 |xm band, and 
compare these data with the new semi-empirical quiet-Sun model from Fontenla 
el al. (2006). The model's solar atmospheric parameters are adjusted to match 
available satellite solar spectral irradiance observations; the model explicitly 
includes the effect of solar activity (both the distribution and radiative effects of 
sunspots, plages, and networks). Model irradiances in the 3.7 pm spectral region 
are expected to have uncertainties of about 2 %; deviations from quite sun values 
are found to be less than 0.5%. 

While a recent study by Trishchenko (2006) examined differences in AVHRR 
and GOES imager 3.7 pm band reflected radiance and brightness temperature 
calculations using four solar irradiance data sets (compiled observational as well 
as synthetic), our current study includes a critical discussion of the original 
observational data as well as the new Fontenla et al. (2006) models, with an 
emphasis on cloud retrieval applications using the MODIS sensor. 

Section 2 gives an overview of the existing 3.7 pm band absolute spectral 
irradiances data sets, in particular, those that have typically been used to convert 
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radiance measurements to solar reflectances for quantitative cloud microphysical 
retrievals. Section 3 summarizes the Fontenla et al. (2006) set of models and 
compares the quiet-Sun spectral irradiance with the historic data sets. The 
resulting differences in the band-averaged solar irradiance for ten AVHRR 
instruments (NOAA-7 through NOAA-18) and both MODIS instruments (Terra 
and Aqua) are given. Relative differences in inferred reflectance using the 
various irradiance data sets are also given for each instrument. We conclude with 
a discussion on the impact of solar irradiance uncertainty on cloud microphysical 
retrievals. 

2. Overview of Historical 3.7 pm Spectral Irradiance Data Sets 

The only observational absolute spectral irradiance data set covering the 
3.7 pm band was obtained during the August 1967 aircraft campaign described 
by Thekaekara et al. (1969, 1974). Instrument and analysis details of the flights off 
the coast of California are extensively discussed in Thekaekara et al. (1969) and 
summarized here. 

Measurements were made from the NASA Convair CV-990 aircraft flying at 
an altitude of 11.5 km. During the six-flight campaign, a variety of total and 
spectral irradiance measurements were obtained from 12 instruments, only two 
of which measured spectral irradiance in the 3.7 pm band: a prism 
monochrometer (operated by Thekaekara et al.) and a Michaelson interferometer. 
To obtain irradiance over the entire solar disk, the monochrometer used a diffuse 
incident mirror with radiometric calibration traceable to a standard lamp; 
corrections were required to account for the additional optical path of sunlight 
through a sapphire window on the aircraft. The interferometer was calibrated to 
an onboard blackbody with corrections for detector and instrument emission. A 
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"weighted average" of irradian.ce from the monochrometer and interferometer 
instruments was reported for tide 3.7 pm band though no details on the 
weighting were given. An overall statement of the derived spectral irradiance 
accuracy was estimated at ±5%. This assessment was not specified as a function 
of wavelength. Spectral irradiance in the 3.7 pm band was reported at 100 ran 
resolution. All irradiances were normalized to the average sun-earth distance (1 
AU). 

Observations were made over a range of solar zenith angles and so 
corrections for slant path (air mass) as well as the atmospheric path absorption 
are critical. In the visible portion of the spectrum, air mass corrections derived 
from the Langley plot method for a single flight day were used to extrapolate 
observations to zero air mass for all flights (with the implicit assumption that the 
above-flight level amount of absorbing gases and aerosol did not significantly 
change during the multi-day observations). Otherwise, it appears that a zero air 
mass correction was obtained from a weighted average of solar zenith angle 
observations and ground-based air mass calculations (details on the weighting 
not described). Atmospheric spectral attenuation calculations are those of 
Elterman and Toolin (1965) for the U.S. Standard Atmosphere. Table 7-4 in 
Elterman and Toolin provides optical properties at 3.5 and 4.0 pm; at both 
wavelengths, the optical thickness between 11 km and TOA is given as zero. 
Molecular absorption in die 3.7 pm band is now known to be important 
(primarily CH 4 and N 2 O in the short and longer wavelength portions of the 
band, respectively, along with water vapor absorption throughout the band). We 
performed our own calculations using the MODTRAN atmospheric 
transmittance code (Berk et al., 2003) for a Midlatitude Summer standard 
atmosphere, an atmospheric path from 11 km to TOA, and an air mass of 2 (air 
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masses over the course of the six CV-990 flights ranged from near-unity to 7). 
Averaged over 50 run bandpasses (estimated grating resolution at these 
wavelengths), the atmospheric path absorptan.ce was 0.024, 0.006, 0.006, 0.040, 
0.021 at 3.6, 3.7, 3.8, 3.9, and 4.0 pm, respectively. Though a critical function of air 
mass, it appears that a significant 2-4% absorption on either side of the 3.7 pm 
atmospheric window was likely even at CV-990 altitudes, and not accounted for 
in the analysis. 

Absolute radiance measurements in the spectral region include those of 
Kondratyev et al. (1965), obtained from Mt. Elbrus in the Caucasus mountains 
(5.6 km altitude) over a two-year period, who reported results for the solar disk 
center at 3.0, 3.6, and 4.0 pm. Their results were converted to irradiance by Pierce 
and Allen (1977), along with the radiance measurements of Koutchmy and 
Peyturaux (1970) from Mt. Louis in the Pyrenees at 3.8 pm and the 31 km altitude 
balloon measurements of Murcray et al. (1964) at 4.0 and 4.1 pm. However, the 
solar limb darkening curve(s) used for the Pierce and Allen calculations was not 
discussed. The derived spectral irradiances for Kondratyev et al. are about 15% 
smaller than Thekaekara et al. at 3.6 pm; the Koutchmy and Peyturaux derived 
irradiance is identical with Thekaekara (at 1 W-m^-pm 4 precision); irradiances 
for Murcray et al. are 14% smaller than Thekaekara et al. at both wavelengths. 
Kondratyev et al. give uncertainties of a few percent in the 3.7 pm band, as does 
Murcray et al. for their 4-5 pm measurements. A summary of these and other 
historical solar spectral observations is given in Pierce and Allen (1977). 

High spectral resolution interferometer observations of solar disk-center 
radiance covering the mid-wave infrared have also been published. These data 
include the ATMOS space shuttle experiment Fourier Transform Spectrometer 
observations (Farmer 1994) and the ground-based Kitt Peak atlas of Livingston 
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and Wallace (1991). However these data were intended for spectroscopic studies 
and are not absolutely calibrated. 

Several compilations of observational and/or modeled solar irradiance 
spectra are in wide use by the remote sensing community. The World Radiation 
Center (WRC) Reference Spectrum (Wehrli 1985) adopted by WMO is a 
compilation from a number of sources. For the present discussion, irradiance in 
the mid-wave infrared portion of the spectrum comes from Smith and Gottlieb 
(1974), which is in turn, a compilation of several sources. All data points from 
Smith and Gottlieb in the 3.7 pm spectral region are from the few Koutchmy and 
Peyturaux (1970) and Murcray et al. (1964) radiance measurements previously 
discussed; again, no information is given regarding the limb darkening curves 
used for the radiance to irradiance conversion in this spectral region. To 
compensate for the sparse data in this spectral region. Smith & Gottlieb provided 
the coefficients of a linear fit to a log-log plot of the available data; this means of 
interpolation is apparently used to provide the 10 nm resolution found in Wehrli. 

A data set by Kurucz (1995) was incorporated into the Moderate Spectral 
Resolution Atmospheric Transmittance (MODTRAN) radiative transfer model 
along with some later modifications (Berk et al. 2003; see file DATA/nezukur.dat in 
the MODTRAN4 distribution). While MODTRAN4 provides options for using 
other irradiance data sets for various parts of the solar spectrum, the Kurucz 
quiet-Sun spectrum is the sole data set covering the 3.7 pm spectral region. From 
the available documentation, the Kurucz spectrum is calculated using a solar 
atmospheric model from Avrett (which we believe is discussed in Fontenla et al. 
1993) that was further modified by Kurucz, and a spectrum computed using lines 
and atomic data compiled by Kurucz. The ASTM E-490 solar spectrum (2000), 
developed for use by the aerospace community, also uses the Kurucz spectra in 
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the 3.7 pm band. 

3. New Spectral Irradiance Calculations 

With only a single absolute observational irradiance dataset available, and 
difficulty and/or lack of information for deriving irradiance uncertainties from 
the solar radiance observations, it is important to include information available 
from solar models in assessing the state of the knowledge for this spectral region. 
Here we give a brief description of calculations using the Fontenla et al. (2006) 
published set of solar models that were further developed by Fontenla et al. 
(2006b) and describe the relevant atomic and molecular data for computing the 
spectrum in the 3.7 pm band. We then compare our results with the historical 
datasets described in Section 2. The newly computed quiet-Sun spectrum from 
Fontenla et al., as implemented in the Solar Radiation Physical Modeling (SRPM) 
system (described below) and covering a spectral range from the UV to 5 pm, is 
planned for inclusion in the upcoming MODTRAN5 release (G. Anderson, 
private communication). 

a. Overview of the model 

The solar irradiance spectrum is calculated from the model atmospheres of 
solar surface features reported by Fontenla el al. (2006). That paper lists a set of 
models of the solar photosphere in which the designations C, E, and F 
correspond to the quiet Sim network (cell center, network lane, and active 
network, respectively). Model C from that dataset is used in this study. Although 
there are differences in the chromospheric parts of these quiet-Sun models, the 
photospheric layers that determine the 3.7 pm irradiance are identical because 
there is no observational evidence for significant effects of quiet-Sun network 
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elements on the irradiance. Solar atmospheric parameters as functions of height 
for the various features observed on the solar disk at mid spatial resolution (~2 
arc sec) are based on earlier work and determined by comparison with available 
ground and space observations. Ground-based measurements of line profiles and 
center-to-limb-variation of the emitted radiance at visible and infrared 
wavelengths provide constraints on the chemical composition and temperature 
gradient in the models. Absolutely calibrated space-based observations of the 
solar spectral irradiance provide an absolute temperature scale. Because the 
observations used for the models were obtained at shorter wavelengths, the 
computed spectral irradiance in the 3.7 pm band relies on the physical model of 
the solar atmosphere derived from those observations and knowledge of 
absorption and emission radiative coefficients in the solar atmosphere (scattering 
is also considered but not very important in this spectral region). 

The emergent radiance spectrum from model C is computed at ten positions 
in the solar disk to determine its angular dependence (equivalent to the center- 
to-limb-variation over the solar disk). Then, the data are combined to produce 
the irradiance spectrum at 1 AU in the manner described by Fontenla et al. (1999, 
2006) and as implemented in the SRPM system. These papers describe model C 
in more detail; its photospheric layers are an evolved version of the Maltby et al. 
(1986) quiet Sun reference model that were modified to better match visible and 
near-infrared observations, especially those from SOLSPEC (Thuillier et al. 2003). 
Further, the current model calculations use recently published elemental 
abundances (as listed in Fontenla et al. 2006) that affect the derivation of the solar 
atmospheric models and the irradiance spectrum. 

The SRPM calculations indicate that the 3.7 pm continuum irradiance 
originates in photospheric layers with gas pressures between 0.57 xlO 5 and 
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1.08x10 s dyne cm -2 (heights from about 10 to 100 km above the T 500 mn=l level, 
and temperatures in the range 5480 to 6310 K), depending on the area of the solar 
disk (higher pressures for the disk center and lower near the limb). This is a 
much broader range than that of formation of the emitted radiance at any given 
position in the solar disk. 

The irradiance spectrum at 3.7 pm consists of a predominant continuum with 
a number of weak CO and OH molecular lines (only a few which are significant). 
Atomic lines of O, N, S, and C are also present but they are not significant for the 
broadband calculations discussed in later sections. The lines we compute are 
evident in Fig. 1 at sub-angstrom resolution, but due to a lack of reliable atomic 
data our calculated spectrum may contain a few less lines than the actual solar 
spectrum. As an example, in the NIST database that is used, there are five lines of 
Mg I around 3.866 pm for which no oscillator strength is provided and therefore 
they are not computed. However, the overall effects of the lines are small and 
only at longer wavelengths (about 4.295 pm, outside the band we consider) do 
CO bands become strong enough to produce a significant reduction in the 
broadband irradiance. 

The main contribution to the opacity in the 3.7 pm band is the free-free 
opacity of the negative ion of Hydrogen. Fontenla, Balasubramanian & Harder 
(2006b) updated older estimates used in Fontenla et al. (2006) to the values 
published by Bell and Berrington (1987) and this produced some changes to the 
infrared intensities computed by SRPM. These current opacity values are thought 
to be accurate to within 2% (John 1988) which, taking into account the 
temperature gradient in the region of continuum formation, would result in 
about a 0.5% irradiance error. Other minor opacity sources at 3.7 pm (in the 
layers that contribute to the irradiance) are the H2 free-free and photoionization 
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of excited states of H with a few percent contribution to the total opacity. 

The SRPM computed irradiance spectrum is 2.2% below Thuillier et al. 
(2003) at their longest measured wavelength (2.30 pm). This suggests that the 
3.7 pm opacity is slightly overestimated in our calculations and thus the 
irradiance is underestimated by perhaps 2%. Such an error is somewhat higher 
than what is estimated from the atomic data uncertainty but well within the 
combined error of the observations (~2% according to Thuillier 2003) and the 
opacity data. 

Fontenla et al. (2006) also include models of features related to solar activity, 
such as sunspots, faculae, plage, and network. Using the computed irradiance 
spectrum for these features coupled with their distribution on the solar disk as 
described by Fontenla and Harder (2005), it was found that the variations in solar 
irradiance in the 3.7 pm spectral region due to solar activity were less than ~0.5% 
of the quiet-Sun irradiance during the last 2 years. The magnitude of these 
variations are confirmed by the variability observed at visible and shorter IR 
wavelengths by the Spectral Irradiance Monitor (SIM) on board the Solar 
Radiation and Climate Experiment (SORCE) satellite (Harder et al. 2005, and 
Rottman et al. 2005). Since the layers producing the 3.7 pm irradiance are the 
same and just a bit shallower than those producing the visible irradiance at the 
continuum around 800 ran, and the observed and modeled irradiance variations 
are close to what is observed at that wavelength, we believe that the model 
estimated variability at 3.7 pm is close to the real variability. These irradiance 
variations are relatively small compared with an uncertainty of about 2% in the 
overall value, which as explained above, is due to uncertainties in the available 
spectral irradiance observations and atomic data. 
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b. Model spectrum and comparisons with other data sets 

Spectral irradiance from Fontenla (2006) is plotted in Fig. 1. Line features are 
evident, especially around 4.05 pm, though they are relatively insignificant 
compared to continuum emission. All subsequent calculations use a reduced 
resolution form of the high-spectral resolution model output shown in the figure 
(convolution of the high resolution data with a 1 run FWHM filter). 

The irradiance curve of Fig. 1 can be accurately fit (mean differences better 
than 0.3 %) over the range 3.40 - 4.15 pm by the following quadratic function: 

F 0 ( A) = 157.91 - 66.34 A + 7.265 A 2 , (4) 

with A in micrometers and Fo in W-m^-pm -1 . Brightness temperature is 
approximately linear across the spectral range (plot not shown) and can be fit to 
the same level of accuracy by 

T b ( A) = 6533.6 - 229.6 A, (5) 


with A in micrometers and Tb in K. 

Plots of the data sets discussed in Section 2 are also shown in Fig. 1. As a 
reminder, other than the Kurucz (1995) spectrum, these historical data are 
observational though all except Thekaedara (1974) are derived from disk center 
radiance measurements. The Thekaekara irradiance is uniformly larger than all 
other datasets except at the shorter wavelengths where values are roughly the 
same as the Fontenla et al. (2006) and Kurucz spectrum. The Fontenla et al. and 
Kurucz spectra generally give irradiances within a couple of percent of each 
other, though the latter model shows more significant line features in the center 
of the 3.7 pm band. There is insufficient documentation on the Kurucz (1995) 
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solar model and atomic and molecular data to compare with those used in our 
calculations. The Whehrli (1985) spectrum is biased high compared to the 
Kontratyev et al. (1965) and Murcray et al. (1964) data points, despite the fact that 
these data were the ultimate source for the Smith and Gottlieb (1974) fits adopted 
by Whehrli; however as noted in the Section 2, the limb darkening curve used by 
Smith and Gotlieb to convert these data to irradiance is not likely the same as the 
curve used by Pierce and Allen (1977) whose tabular values are shown in Fig. 1. 

c. Instrument-specific sensitivity to the irradiance spectrum 

Uncertainty in inferring the bi-directional reflectance from a 3.7 pm channel 
observation is dependent on the uncertainty in the instrument spectral response 
function, the solar irradiance spectrum, and the instrument's radiance calibration 
uncertainty (Eq. 1). 

The relative spectral response functions for the 3.7 pm channel for four of the 
ten AVHRR instruments considered in this study, and the two MODIS 
instruments, are shown in Fig. 2. The Fontenla and Thekaekara spectral 
irradiance data (from Fig. 1) are also shown. The calculated bandpass and central 
wavelength for both MODIS and all ten AVHRR instruments is given in Table 1. 
The two MODIS response functions are very similar while there is wide variation 
in the AVHRR channel characteristics (differences in central wavelength of over 
90 nm between NOAA-12 and NOAA-15). All response functions are based on 
pre-launch measurements. While both instruments are filter radiometers, MODIS 
interference filters are made with a more stable ion-assisted deposition process 
(Montgomery et al. 2000). Further, MODIS has an onboard calibration source to 
monitor in-orbit changes in the channel center for spectral channels from the 
visible up to 2.1pm (the Spectral Radiometric Calibration Assembly, 
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Montgomery et al.) and no significant spectral shifts have been observed in these 
channels for either MODIS instrument (e.g., less than 0.2 ran for shortwave 
infrared channels, Xiong et al. 2005a). While 3.7 pm spectral stability cannot be 
monitored in-orbit, performance similar to the shorter wavelength channels is 
expected. 

The differences in filter location and bandpass seen in Fig. 2 result in 
substantial differences in the calculated average solar spectral irradiance over the 
channel (Eq. 2); such calculations for each instrument are given in Table 2 for the 
Fontenla (2006), Thekaekara (1974), and Kurucz (1995) spectra. The differences in 
channel-averaged irradiances relative to the Fontenla spectrum are given by the 
first pair of data columns in Table 3. As an example, for the MODIS Terra 
response function, use of the Thekaekara data results in an averaged spectral 
irradiance 3.85% larger than that obtained from Fontenla and 5.45% larger than 
use of the Kurucz spectrum. Similar values are found in the second pair of 
columns in Table 3 where the irradiance differences have been given in terms of 
solar reflectance (Eqs. 1, 2). Note that the difference in the average irradiance 
between the two MODIS instruments is less that 0.5%, regardless of the 
irradiance spectrum used, while AVHRR differences can be in excess of 6% 
(NOAA-12 vs. NOAA-15). 

For context, the calibration of MODIS solar bands used for cloud retrieval 
algorithms (visible through the 2.1 pm band) is derived from an onboard solar 
diffuser; including transfer of the diffuser calibration to the onboard diffuser 
observation, uncertainty is less than 2% for most solar bands (Xiong et al. 2005b). 
Instriunent-specific uncertainties related to the 3.7 pm characterization of 
onboard blackbody calibration sources and their transfer need to be added 
appropriately to the numbers in Table 3 (e.g., addition of variances assuming 
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zero correlation) to provide an overall reflectance measurement uncertainty. The 
total RMS onboard radiance calibration uncertainty for the MODIS 3.7 pm 
channel (referred to as band number 20) has been evaluated to be less than 0.7 % 
over a wide rand of radiances for both Terra and Aqua instruments (Xiong et al. 
2005b). Therefore, for MODIS, the solar irradiance uncertainties of Table 3 
dominate the radiometric and spectral uncertainties and thus also dominate the 
reflectance calculation uncertainty. Issues related to AVHRR thermal channel 
radiometric uncertainty and stability are discussed by Trishchenko (2002). 

d. Cloud retrieval sensitivity 

The significance of the reflectance differences in Table 3 for cloud retrieval 
applications was empirically investigated using a Terra MODIS data granule of 
coastal Chile and Peru that has been previously discussed (Platnick et al. 2003). 
The granule (5 minutes of data) includes a variety of cloud types, including 
marine boundary layer water clouds, cirrus over land and ocean, and convective 
ice and water clouds over the Amazon basin. Retrievals were obtained from the 
operational MODIS Terra cloud product (MOD06) Collection 5 processing 
stream code, that includes separate effective retrievals derived from 3.7 nm 
MODIS observations (as well as size retrievals derived from 1.6 and 2.1 jxm 
observations, separately and in combination). Histograms of liquid water and ice 
cloud effective radius retrievals for the data granule are shown in Fig. 3 using 
two different values for the 3.7 pin MODIS channel-averaged spectral irradiance 
(a 5% difference). The 11.34 W-nr 2 ~pm 4 value corresponds to the average of the 
MODIS Terra and Aqua instrument irradiances using the Thekaekara (1974) data 
set. The smaller irradiance (10.77 W-nr 2 -|xm 4 ) is comparable to the instrument 
averages using either Fontenla (2006) or Kurucz (1995). For liquid water clouds. 
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the difference in the mean value of the retrieved effective radius over the data 
granule is 0.50 [im (»4.5 %); for ice clouds, the difference is 0.7 1 pm (« 6.5 %). 

To place these differences in context, the 3.7 fim effective radius retrieval 
uncertainty for these MODIS retrievals was calculated assuming the spectral 
irradiance is known exactly, but including the following set of error sources 
contributing to the reflected signal: instrument calibration and forward model 
uncertainty (5%), surface spectral albedo uncertainty (15 %), above-cloud column 
water vapor amount derived from model analysis and MODIS cloud-top 
pressure retrievals for use in spectral atmospheric corrections (20%). Each error 
source is assumed uncorrelated. Note that the uncertainty analysis does not 
include errors in the estimation of cloud-top temperature, surface temperature, 
and atmospheric state that are required to account for the 3.7 pm emission 
component. The retrieval uncertainty is a function of effective radius and cloud 
optical thickness, as well as surface type and atmospheric state, and as such, is 
difficult to summarize succinctly (Platnick et al. 2004). However, for comparison 
with the granule-averaged means of Fig. 3, mean 3.7 pm effective radius relative 
uncertainties are 5.2% and 6.4% for liquid water and ice clouds, respectively. 
Therefore, in this example, the uncertainty in the retrieved effective radius due to 
a 5% irradiance uncertainty is of the same order as other fundamental reflectance 
error sources. 

4. Discussion 

Despite the widespread presence of a 3.7pm atmospheric window channel 
on present-day and planned satellite imagers, the absolute top-of-atmosphere 
solar spectral irradiance in the band, which is required for daytime quantitative 
use, is limited to a single aircraft field campaign (Thekaekara et al. 1969, 
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Thekaekara 1974). This sole observational irradiance data set was discussed, 
along with a survey of published radiance, observations, and compared with 
historical synthetic data sets and the new solar models of Fontenla et al. (2006, 
2006b). Channel-averaged spectral irradiances using the various data sets were 
calculated for laboratory-derived AVHRR and MODIS instrument spectral 
response functions. Depending on the particular instrument, the observations 
reported by Thekaekara (1974) give averaged irradiances from about +3.2% to 
+4.1 % greater than the quiet-Sun Fontenla et al. spectrum; the Kurucz spectra 
gives channel-averaged irradiances about -1.2% to -1.5% smaller than Fontenla et 
al. The new Fontenla et al. irradiance spectrum is expected to have a 2% 
uncertainty in the 3.7 pm region. Accurate analytic fits to the Fontenla et al. 
spectral irradiance and brightness temperature across typical 3.7 and 3.9 pm 
channels have been provided. The quiet-Sun irradiance spectrum from Fontenla 
et al. used in this study (as implemented in the Solar Radiation Physical 
Modeling (SRPM) system) is planned for inclusion in the upcoming MODTRAN5 
release. 

The consequences on cloud microphysical retrievals from uncertainty in the 
3.7 pm band solar irradiance was examined for a MODIS data granule consisting 
of a variety of liquid and ice clouds over both ocean and land. Using a 
representative channel-averaged solar irradiance uncertainty of 5%, based on 
this study's comparisons of observational and modeled irradiances, averaged 
retrieved effective radius uncertainties averaged for all pixels in the data granule 
were 4.5% and 6.5% for liquid water and ice clouds, respectively. This is 
comparable to retrieval uncertainties due to other error sources that 
fundamentally impact the inference of cloud-top reflectance in the band. 

While the relative agreement between the Fontenla et al. (2006) spectrum 
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and Thekaekara (1974) in the 3.7 um band is reassuring, new absolute spectral 
irradiance observations in the band are sorely needed to augment the sole 
available observations reported by Thekaekara. 
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Table 1. The spectral response function Full- Width-Half-Maximum (FWHM) 
bandpass and central wavelength (defined as the average of the two FWHM wavelengths) 
for AVHRR and MODIS 3.7 pm channels used in this study. 


Instrument/ 

Platform 

Central Wavelength 
(pm) 

Bandpass 

(pm) 

AVHRR N7 

3.741 

0.411 

AVHRR N9 

3.733 

0.410 

AVHRR N10 

3.763 

0.375 

AVHRR Nil 

3.750 

0.400 

AVHRR N12 

3.800 

0.400 

AVHRR N14 

3.788 

0.425 

AVHRR N15 

3.708 

0.322 

AVHRR N16 

3.724 

0.326 

AVHRR N17 

3.761 

0.385 

AVHRR N18 

3.768 

0.385 

MODIS Terra 

3.792 

0.190 

MODIS Aqua 

3.785 

0.189 


i 
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Table 2. Channel-averaged solar spectral irradiance for AVHRR and MODIS 3.7 pm 
instrument channels, for three spectral irradiance data sets. 


Instrument/ 

Platform 

Band-Averaged Solar Spectral Irradiance (W-m^-pm' 1 ) 
Solar Irradiance Dataset 

Fontenla et al. 
(2006) 

Thekaekara 

(1974) 

Kurucz 

(1995) 

AVHRR N7 

11.573 

11.957 

11.429 

AVHRR N9 

11.671 

12.043 

11.528 

AVHRR N10 

11.360 

11.771 

11.216 

AVHRR Ni l 

11.543 

11.930 

11.400 

AVHRR N12 

11.020 

11.470 

10.879 

AVHRR N14 

11.138 

11.573 

10.997 

AVHRR N15 

11.729 

12.101 

11.577 

AVHRR N16 

11.467 

11.867 

11.317 

AVHRR N17 

11.327 

11.741 

11.184 

AVHRR N1 8 

11.200 

11.627 

11.059 

MODIS Terra 

10.885 

11.304 

10.720 

MODIS Aqua 

10.974 

11.386 

10.807 
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Table 3. Difference in the 3.7 pm channel-averaged solar spectral irradiance and 
inferred solar reflectance relative to the Fontenla et al. (2006) data set, as derived from 
Table 2. 


Instrument/ 

Platform 

Relative Difference (%) in Solar 
Irradiance vs. Fontenla 

Irradiance Dataset 

Relative Difference (%) in Inferred 
Solar Reflectance vs. Fontenla 

Irradiance Dataset 


Thekaekara 

(1974) 

Kurucz 

(1995) 

Thekaekara 

(1974) 

Kurucz 

(1995) 

AVHRRN7 

3.32 

-1.24 

-3.21 

1.26 

AVHRRN9 

3.19 

-1.23 

-3.09 

1.24 

AVHRRN10 

3.62 

-1.27 

-3,49 

1.28 

AVHRRN11 

3.35 

-1.24 

-3.24 

1.25 

AVHRRN12 

4.08 

-1.28 

-3.92 

1.30 

AVHRRN14 

3.91 

- 1.27 

- 3.76 

1.28 

AVHRRN15 

3.17 

-1.30 

- 3.07 

1.31 

AVHRRN16 

3.49 

-1.31 

-3.37 

1.33 

AVHRRN17 

3.65 

-1.26 

-3.53 

1.28 

AVHRRN18 

3.81 

-1.26 

- 3.67 

1.27 

MODIS Terra 

3.85 

-1.52 

-3.71 

1.54 

MODIS Aqua 

3.75 

-1.52 

-3.62 

1.55 
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Fig. 1 . Plot of selected solar spectral irradiance datasets discussed in the text for the 
3.7 pm band spectral region. Data points for Kondratyev et al. and Murcray et al. were 
taken from Pierce and Allen (1977). 
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2. 3.7 pm channel relative spectral response functions for the two MODIS instruments (NASA Terra and Aqua spacecraf 
et of the AVHRR instruments (NOAA polar orbiting platforms 9, 14, 15, and 16) used for channel-averaged irra 
itions. The solar spectral irradiance from Fontenla et al. (2006) and Thekaekara (1974), taken from Fig. 1, is also shown. 
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Fig. 3. Frequency histogram for liquid water and ice cloud retrievals from a MODIS 
Terra data granule (see text for details) using two different values for the channel- 
averaged spectral iiradiance. The mean effective radius for the retrievals is given in 
parenthesis. 





